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ABSTRACT 

Plane Poiseuille flow shows turbulence at a Reynolds 
number that is lower than the critical one for the onset of 
Tollmien-Schlichting waves. The transition to turbulence 
follows the same route as the by-pass transition in boundary 
layers, i.e. finite amplitude perturbations are required and 
the flow is dominated by downstream vortices and streaks in 
the transitional regime. In order to relate the phenomenol¬ 
ogy in plane Poiseuille flow to our previous studies of plane 
Couette flow ( [Kreilos & Eckhardt] |2012| ), we study a sym¬ 
metric subspace of plane Poiseuille flow in which the bi¬ 
furcation cascade stands out clearly. By tracing the edge 
state, which in this system is a travelling wave, and its bi¬ 
furcations, we can trace the formation of a chaotic attractor, 
the interior crisis that increase the phase space volume af¬ 
fected by the flow, and the ultimate transition into a chaotic 
saddle in a crisis bifurcation. After the boundary crisis we 
can observe transient chaos with exponentially distributed 
lifetimes. 


Introduction 

The laminar state of plane Poiseuille flow (PPF) be¬ 
comes linearly unstable at a Reynolds number of 5772.22 
( |Orszag[ \\91l) . Above this critical Reynolds number 
any small disturbance lead to the formation of Tollmien- 
Schlichting (TS) waves. As time progresses, these waves 
grow in amplitude until they finally undergo a secondary in¬ 
stability, break up, and a reach a turbulent state. Below the 
critical Reynolds number, just as for pipe and plane Couette 
flow, turbulence can be observed down to Reynolds num¬ 
bers between 600 and 1000 ( [Carlson et al.\ 11982) |Tucker^] 
|man et al\ |2014[ |Xiong et ~cd\ 2015| ) although the laminar 
state is linearly stable. Moreover, while the transition mech¬ 
anism via TS-waves is subcritical, the subcritical branch 
does not reach sufficiently low to explain the transition in 
that range. Therefore, the transition mechanism has to be 
different from that associated with TS-waves. 


Much insight into the process can be gained by study¬ 
ing the critical amplitudes that have to be surpassed in order 
to trigger transition. They depend on Reynolds number and 
decrease with increasing Reynolds number ( |Lemoult et al. \ 
2012). More importantly, the boundary between initial con¬ 
ditions that decay and those that reach turbulence is formed 
by states that are on the stable manifold of an intermediate 
state, the edge state of the system ( |Skufca et a/.[|2006) . The 
different states and their variation with Reynolds number 
are shown in figure [l] 


The edge state also provides access to other coherent 
structures relevant for the dynamics. Tracing it to lower 
Reynolds number one reaches the point where the state ap¬ 
pears in a saddle-node like bifurcation. Following the up¬ 
per branch back to higher Reynolds numbers then gives 
states that characterize the chaotic dynamics, i.e. the transi¬ 
tion from simple exact coherent structures to more compli¬ 
cated ones, very often in standard bifurcation such as Hopf- 
bifurcations and period doubling cascades. This connection 
between exact coherent structures and the temporal prop¬ 
erties of the chaotic dynamics was studied in plane Cou- 
ette ([Kreilos & Eckhardt[|2012] > and pipe flow ( [Avila et al\ 
|2013| ). For the case of plane Couette flow it was found that 
the upper branch of the Nagata-Busse-Clever-solution, the 
first three-dimensional equilibrium solution identified for 
plane Couette flow ( |Nagata[[l990[|Clever & Busse||1992| ), 
undergoes a series of bifurcations leading to a chaotic at¬ 
tractor. A boundary crisis bifurcation ( [Lai & T dT| |2011| ) 
turns the attractor into a chaotic saddle. Following the 
boundary crisis, the lifetimes are exponentially distributes, 
as observed in experiments ( Bottin & Chate 1998) and sim¬ 
ulations ( Schnei der et al\ |2010| >. In pipe flow the same 



Figure 1. Bifurcation diagram of plane Poiseuille flow. 
The laminar state becomes unstable to TS-waves in a sub¬ 
critical bifurcation at Re = 5772.22. The TS-waves ex¬ 
ist down to Re ~ 2600 for the optimum wavelength (light 
blue). At low Reynolds numbers the edge state appears and 
bypass-transition is possible for finite amplitude perturba¬ 
tions. The critical Reynolds number for the appearance of 
exact solutions is typically lower than the Reynolds number 
where turbulence is observed in experiments, because near 
the critical point only a small part of state space is affected. 
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mechanism exist. Here, the relevant state that is the start¬ 
ing point of the cascade is a streamwise localized periodic 
orbit which can be found using the technique of edge track¬ 
ing (|Skufca et al.\ |2006J in a symmetric subspace (|Avila| 
|efaz'||2013). ~ 


We recently found the same mechanism for the genera¬ 
tion of subcritical turbulence in PPF jZammert & Eckhardt[ 
|2015f . For a symmetric subspace of PPF we will present 
the edge state of the system and discuss the formation of 
the chaotic saddle. Furthermore, we will show that attrac¬ 
tors and saddles with quite different spatial wavelength exist 
in the state space and discuss possible implications of this 
observation. 

For our numerical simulation of plane Poiseuille flow 
we use the Open-Source Channelflow-code developed by 
J.F.Gibson ( |Gibson[[2012j ). The code uses periodic bound¬ 
ary conditions in streamwise and spanwise direction to¬ 
gether with no-slip boundary conditions at both walls. A 
constant mass-flux is enforced in all our calculation. The 
definition of the Reynolds number Re = U^d V in this work 
uses half the distance between the plates d and the center- 
line velocity Uq of the laminar profile. 


The edge state in a symmetric subspace 

We study PPF with two imposed symmetries in a small 
periodic domain with a streamwise and spanwise wave¬ 
lengths of 2tz and 7T, respectively. The used numerical res¬ 
olution is N x x N y x N z = 48 x 65 x 48. The imposed sym¬ 
metries are a reflection at the center-plane, 


Sy : [u,v,w]{x,y,z) = [u,-v,w](x,-y,z), (1) 


and a reflections in spanwise direction, 


s z : [u,v,w\{x,y,z) = [«, v, -w](x,y, -z). (2) 


In this symmetric subspace the dynamics is simpler and ex¬ 
act coherent structures within this subspace are stabilized, 
since all unstable direction pointing out of the subspace are 
removed. 

We recently showed ( [Zammert & Eckhardt]|2015| ) that 
in this symmetry restricted subspace the edge state, the at¬ 
tracting state on the boundary between laminar and turbu¬ 
lent initial conditions, is a travelling wave, henceforth re¬ 
ferred to as TWE yz . This travelling wave can be identified 
using edge tracking (Skuf ca et~ai\ [2006). The technique 
uses a bisection in amplitude to find trajectories on both 
sides of the laminar-turbulent boundary as approximation 
to a trajectory on the laminar turbulent boundary. To stay 
close to the laminar-turbulent boundary from time to time 
refinements by additional bisections are necessary. Using 
such refinement steps, edge tracking allows to follow tra¬ 
jectories on the laminar turbulent boundary for an arbitrary 
long time and thus to find attracting states in the boundary. 

In figure [2] the edge tracking used to find TWE yz is 
shown. The edge trajectory, approximated by trajectories 
on the laminar (blue) and the turbulent (red) side of the edge 
quickly reaches a state of constant energy 
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Figure 2. Edge trajectory at Re = 2000 in the symmetric 
subspace. For each refinement the red and blue lines are 
used for the trajectories on the turbulent and laminar side, 
respectively. Black circles mark the positions of refinement 
steps, additional bisections, necessary to stay close to the 
trajectory on the edge. 


This indicates that the attracting state is a relative equi¬ 
librium. Using an initial condition on the edge trajectory 
as initial guess, a Newton-search for a relative equilibrium 
converges to the travelling wave TWE yz 

The wave can be tracked in Reynolds number using 
a numerical continuation technique what allows to find its 
turning point and the upper branch can. Visualizations of 
the lower and upper branch of the traveling wave are shown 
in figure [3] The visualizations using iso-contours of the 
streamwise velocity (deviation form laminar) and the Q- 
vortex criterion (Jeong & Hussain 1995) show that the state 
consists of streamwise streaks and vortices. On the lower 
branch the vortices are located near the center-line while 
on the upper branch the strongest vortices can be found 
closer to the walls. Furthermore, iso-contours of the Q- 
vortex criterion suggest that on the upper branch the vor¬ 
tices have a hairpin-like structure because the tips of two 
counter-rotating vortex tubes connect. 


The formation of a chaotic saddle 

A series of bifurcations that start at the upper branch of 
the state at Re = 707 creates a chaotic attractor. The series 
starts with the bifurcation of a relative periodic orbit and 
involves different kinds of bifurcations, e.g. period dou¬ 
blings and interior crisis bifurcations (Grebogi et al. 1982, 
Lai & Tel 2011). Ultimately, the attractor is converted into 
a chaotic saddle by a boundary crisis bifurcation. Sampling 
the minima and maxima of the energy along trajectories on 
the attractor it is possible to visualize the cascade in a bi¬ 
furcation diagram. In figure [4] the diagram is shown for 
Reynolds numbers close to boundary crisis. 

To investigate the influence of the boundary crisis on 
the state space structure, two-dimensional slices of the state 
space are regarded. Using two parameters, a and A, and the 
flow fields of the lower (u\) and upper (U 2 ) branch of TWE yz 
a slice of the state space can be defined by: 


u(a,A) =A 


(1 — a)u\ + au 2 
\/E(( 1 — a)u\ + 0 LU 2 ) 


(4) 


(3) 
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The parameter a allows for interpolation between the lower 
and upper branch and A is the amplitude of the flow field. 
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Figure 4. Bifurcation diagram of the attractor in a compu¬ 
tational domains with L x = 2k and L z = Itt. The states are 
visualized by plotting minima and maxima of E(u) along 
trajectories. Continues lines then indicate states with a sim¬ 
ple periodicity, splittings (as near 784) an increase in period 
and complexity, and clouds an irregular temporal variation. 
The red line markes the Reynolds number where the fluc¬ 
tuations in the state meet the lower branch state and the at¬ 
tractor breaks up into a transient chaotic saddle. 


Figure 3. Visualizations of the lower (a) and upper branch 
(b) of the travelling wave TWE yz at Re = 850. The red and 
blue surfaces show iso-surfaces of u = ±0.15 for the lower 
and u = ±0.30 for the upper branch. The yellow surfaces 
show iso-surfaces of the Q-vortex criterion. The iso-values 
are 0.008 and 0.08, respectively. On the back side of the box 
the streamwise velocity (deviation from laminar) is color 
coded. The vortices in the lower branch are fairly straight 
with weak modulations in spanwise direction. The upper 
branch has a more complex structure. 

Slices for two different Reynolds number are shown 
in figure [5] In these slices the initial conditions are color- 
coded according to the time needed to approach the lami¬ 
nar state. (In practical we use the time the initial condition 
needs to fall below an certain very low value in energy.) 
Below the boundary crisis (figure [5^) initial condition are 
either captured by the attractor (dark red) or return imme¬ 
diately to the laminar state (blue). After the crisis (figure 
[5» the lifetime landscape is fractal as also observed in ear¬ 
lier investigations of plane Couette flow b y|Schmiegel &| 
|Eckhardt| ( |1997| > and |Kreilos & EckhardtlffiT 2]). The distri- 
bution of the lifetimes is exponential with a characteristic 
timescale that decreases with increasing Reynolds number 
and diverges with approaching the Reynolds number of the 
crisis. 

It should be remarked that the wave TWE yz does not 
connect to the Tollmien-Schlichting waves bifurcating from 
the laminar state. Thus, the occurrence of subcritical tur¬ 
bulence in PPF is independent of the linear instability and 
completely analogous to the case of plane Couette or pipe 
flow. 


Attractors and saddles with different wave¬ 
length 

The attractor and the saddle discussed in the previous 
section have a streamwise and spanwise wavelength of 2k 
and 17T, respectively. The Reynolds numbers where the at¬ 


tractor appears and where it is converted into a chaotic sad¬ 
dle depend, of course, on these spatial wavelengths. How¬ 
ever, the phenomenon itself is quite robust to changes of 
these wavelengths and can be observed over a larger range 
of parameter values. 

For streamwise wavelengths of l.5n and 3k we study 
the attractors more detailed. In figure [6] both attractors are 
shown together with the one for L x = 2k. The details of the 
bifurcation diagrams are strongly different but each attrac¬ 
tor finally disappears in a boundary crisis bifurcation. The 
Reynolds numbers of the boundary crisis for 1.5 K and 3 k 
are 781 and 854, respectively. 

Although, the dynamics of the attractor created in the 
bifurcation cascade for the case L x = 2k is chaotic the com¬ 
plexity of the flow structure is rather simples due to the sym¬ 
metry constrain and the small size of the periodic domain. 
The merging of attractors with different spatial wavelength 
(e.g in a merging crisis) is a possible scenario that could 
add spatial complexity to the trajectories. In a computa¬ 
tional domain with a streamwise length of 6k and a width 
of 2k all three attractors or saddles, respectively, are present 
and mergings could actually exist. 


Conclusion 

Subcritical turbulence in plane Poiseuille flow is gen¬ 
erated in a bifurcation cascade starting at a the simple trav¬ 
elling wave TWE yz . Although, the travelling wave and 
the cascade were identified and studied using a symmetry 
restricted system, the are also present in the full system 
but in an unstable subspace. TWE yz does not connect to 
the TS-wave or the laminar flow. Thus, the state and the 
boundary formed by the stable manifold of its lower branch, 
which separates two parts of the state space, can also exist 
above Re cr i t . This might explain why bypass-transition and 
TS-transition coexist in the supercritical range. For a de¬ 
tailed understanding of the coexistence of both transitions 
the symmetry restriction must be dropped in further stud¬ 
ies, because the TS-wave do not exist in the studied mirror- 
symmetric subspace. Since coexisting transition scenarios 
are a phenomenon that also exists in other flows with a base 
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Figure 5. State space slices for different Reynolds num¬ 
bers. The initial conditions are color coded according to 
their lifetime, the time they need to fall below 10 -4 in 
energy. Initial condition that do not fall below this value 
within the observation time of 10 4 time units are shown with 
this value as lifetime. For a) and b) the Reynolds numbers 
are 780 and 790, respectively. The positions of the upper 
and lower branch of TWE yz are marked with black dots. 
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Figure 6. Bifurcation diagrams for states with spatial 
wavelength of 3ti (red), 2n (black) and l.5n (blue). A solid 
line is used for the travelling wave TWE yz . The time depen¬ 
dent states are visualized by plotting minima and maxima 
of E(u) along trajectories. 


flow that has a linear instability, such as the Blasius- or 
the asymptotic suction boundary layer, the results found for 
PPF can also improve our understanding of these flows. 

At low Reynolds numbers PPF (Tuckerman et al. 


2014]) and other su bcritical flows (e.g. Prigent et al. \ |2003 


Moxey & Barkley[|2010[|Duguet et al. 2010] ) show a com- 

plex spatio-temporal behavior . The described existence of 
saddles with different spatial wavelengths and the possibil¬ 
ity of merging might provide a possible to understand this 
phenomenon from a dynamical systems point of view. 


This work was supported in part by the DFG within 
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